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ABSTRACT 

I propose that solar coronal heating is a self-regulating process that keeps the coronal plasma 
roughly marginally coUisionless. The self-regulating mechanism is based on the interplay of two 
effects. First, plasma density controls coronal energy release via the transition between the slow 
collisional Sweet-Parker regime and the fast coUisionless reconnection regime. This transition takes 
place when the Sweet-Parker layer becomes thinner than the characteristic coUisionless reconnection 
scale. I present a simple criterion for this transition in terms of the upstream plasma density (ne), 
the reconnecting (Bq) and guide {Bz) magnetic field components, and the global length (L) of the 
reconnection layer: L < 6 ■ 10^ cm (ne/10^°cm~'^)^^ (i?o/30 G)'' (Bq/Sz)^. Next, coronal energy 
release by reconnection raises the ambient plasma density via chromospheric evaporation and this, 
in turn, temporarily inhibits subsequent reconnection involving the newly-reconnected loops. Over 
time, however, radiative cooling gradually lowers the density again below the critical value and fast 
reconnection again becomes possible. As a result, the density is highly inhomogeneous and intermittent 
but, statistically, does not deviate strongly from the critical value which is comparable with the 
observed coronal density. Thus, in the long run, the coronal heating process can be represented by 
repeating cycles that consist of fast reconnection events (i.e., nanoflares), followed by rapid evaporation 
episodes, followed by relatively long periods ( < 1 hr) during which magnetic stresses build up and 
simultaneously the plasma cools down and precipitates. 

Subject headings: magnetic fields — MHD — Sun: corona — Sun: magnetic fields — Sun: fiares 



1. INTRODUCTION 

In this paper I address some aspects of solar coronal 
heating (see Aschwanden et al. 2001b and Klimchuk 2006 
for recent reviews) in the context of Parker's nanoflare 
model (Parker 1972, 1983, 1988). A more concise version 
of this work is presented in Uzdensky (2006, 2007). 

Since the main heating process in the nanofiare model 
is magnetic reconnection, I first discuss what we have 
learned about reconnection in the past 20 years or so 
(see §[!]). I argue that, even though we still do not have 
a complete picture of reconnection, there now appears to 
be some consensus in the magnetic reconnection commu- 
nity about some of its fundamental aspects. One of the 
main goals of this paper is to use this emerging knowl- 
edge to shed some new light on the old coronal heating 
problem. In this paper I purposefully adopt a rather 
conservative approach: I invoke only those very few re- 
sults that seem to be relatively firmly established and try 
not to rely on those details of reconnection physics that 
are still under vigorous debate. Specifically, there is now 
strong evidence coming from numerous numerical simu- 
lations and some laborious laboratory experiments that 
there are two main modes of reconnection: slow Sweet- 
Parker reconnection taking place in collisional plasmas 
where classical resistive MHD applies; and fast Petschek- 
like reconnection in coUisionless plasmas (in § 12. ip . The 
transition between these two regimes seems to be rather 
sharp. An approximate condition for this transition 
can be formulated as a relationship between the global 
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length L of the reconnecting system and the electron col- 
hsional mean-free path inside the layer, Ae,mfp (see i) l2.2p . 
Furthermore, this condition can be cast in terms of the 
plasma density n, the reconnecting magnetic field Bq, 
and L; namely, one can define a critical density nc{L, Bq) 
below which reconnection switches to the fast regime. 
In the strong-guide field case, B^ 3> Bq, the condition 
is modified; namely, the critical density for transition 
to fast reconnection is suppressed by a factor of order 
{Bz > -Bo) (see § [O]). In § [23] I address various 
alternative ideas and caveats that may complicate the 
above simple picture. 

I then apply these results to the active solar corona 
(in §[3|), viewed withing the nanoflare model. My main 
point here is that the corona should be regarded as a self- 
regulating machine keeping itself (in a statistical sense) 
around marginal coUisionality. This conclusion comes 
from the interplay between the way the plasma density 
controls reconnection via the above coUisionless recon- 
nection transition, and the way the coronal magnetic en- 
ergy release due to reconnection in turn controls the am- 
bient gas density via chromospheric evaporation. The 
coronal heating process is then highly intermittent and 
inhomogeneous; it can be thought of a sequence of char- 
acteristic energy-circulation cycles that occur simulta- 
neously in a broad range of spatial, temporal, and en- 
ergy scales. Each such elementary cycle consists of sev- 
eral phases: (1) a fast reconnection event (a nanoflare) 
causing (2) an evaporation episode filling the loop with 
hot dense plasma, followed by (3) a longer period during 
which the magnetic stresses build up and the plasma den- 
sity goes down due to slow radiative cooling (and thermal 
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conduction) . It is the coUisionless reconnection condition 
that makes this scenario (in particular, the last phase) 
possible. Indeed, an important point here is that, even if 
a current sheet is formed, it will stay in the slow Sweet- 
Parker state if the density is large enough. This will 
continue for some time, until the reconnecting magnetic 
field becomes strong enough and/or the ambient plasma 
density becomes low enough for the system to transi- 
tion to the fast coUisionless reconnection regime. This 
enables a non-trivial amount of free magnetic energy to 
be accumulated before a sudden release. This energy 
amount, along with the characteristic time-scale between 
nanoflares and the characteristic coronal density, is de- 
termined, statistically, by the balance between the rate 
at which current sheets are created and amplified by the 
photospheric footpoint motions and the efficiency of ra- 
diative cooling. The coUisionless reconnection condition 
plays a key role in this process as a mediator and ulti- 
mately governs the statistical distribution of the coronal 
reconnection events (flares). 

In § HI I discuss some of the open questions that, in 
my view, need to (and can) be addressed in the near 
future, in order to see whether the physical picture put 
forward in this paper is correct and what modifications 
should be made to improve it. This section is mostly tar- 
geted towards researchers doing numerical simulations of 
reconnection and also to experimentalists and observers. 
Finally, I present my conclusions in § [5j 

I also would like to make a clarifying remark about 
the terms "coUisional" and "coUisionless" that will be 
used many times throughout this paper. Often, by 
"coUisionless" one means a regime in which fie > ly^^i, 
i.e., when the electron gyro-frequency is higher than the 
electron-ion collision rate. Equivalently, this means that 
Ae,mfp > Pe- In this paper, on the other hand, the 
terms "coUisional" and "coUisionless" will involve com- 
paring Ae,mfp not with the gyro-radius, but instead with 
the system size L, or, more precisely, with L divided by 
the square root of the ion/electron mass ratio. As will 
be explained in §§ 12.2112.31 [specifically, see cq. (U) and 
the discussion at the end of § 12. 3j , this choice of termi- 
nology turns out to be more relevant to the reconnection 
problem. 

2. FAST COLLISIONLESS RECONNECTION CONDITION 

2.1. The Existence of Two Reconnection Regimes 

Magnetic reconnection research started 50 years 
ago with the Sweet-Parker (SP) theory (Sweet 1958; 
Parker 1957) for solar flares. As was realized immediately 
at that time, this relatively simple and elegant theory 
could not reproduce the very short 10^ sec) reconnec- 
tion timescales required by flare observations; instead, 
it predicted reconnection times of order a few weeks or 
months. After several years, however, Petschek (1964) 
proposed a modification to the classic Sweet-Parker the- 
ory that apparently resulted in a much higher reconnec- 
tion rate, thus eliminating the main contradiction with 
observations. He realized that the main bottleneck in the 
Sweet-Parker reconnection model is the need to have a 
reconnection layer that is both thin enough for the resis- 
tivity to be important and at the same time thick enough 
for the plasma to be able to flow out. The result of this 
compromise is the famous Sweet-Parker scaling for the 
thickness (5sp of the reconnection layer and for the recon- 



nection velocity v^cc'- 

where L is the global length of the reconnection layer, 
of order the system size, "q is the magnetic diffusivity, 
and Va is the Alfven speed corresponding to the recon- 
necting (in-plane) magnetic field component. Since the 
Lundquist number in the solar corona is usually very 
large (10^^ or greater), the layer is very thin and the 
resulting reconnection rate is very small. Furthermore, 
Petschek (1964) proposed that this difficulty can be cir- 
cumvented if the reconnection region has a certain special 
structure: the famous Petschek conflguration, with four 
slow-mode shocks attached to a small Sweet-Parker cen- 
tral diffusion region. As he showed, this structure gave 
an additional geometric factor that could lead to a much 
faster reconnection rate. 

I would like to emphasize the special importance of 
Petschek's idea for reconnection in astrophysical sys- 
tems, including the solar corona. As I will discuss below, 
there are several local physical effects (e.g.. Hall effect 
or anomalous resistivity) that can indeed broaden the 
reconnection layer and prevent it from collapsing down 
to the very small Sweet-Parker thickness ^sf- However, 
each of these processes comes with its own microscopic 
physical scale, e.g., the ion gyro-radius pi or the ion 
coUisionless skin-depth di = c/tUpi. These scales are 
determined by the local values of the basic plasma pa- 
rameters inside the reconnection layer, such as magnetic 
fleld, density, or temperature. The important point is 
that they know nothing about the overall global system 
size L. Astronomical systems are often astronomically 
large, however. That is, L is typically much larger than 
any microscopic physical scale S and than the Sweet- 
Parker layer thickness (which is a hybrid length-scale). 
Therefore, any simple Sweet-Parker-like analysis would 
give a reconnection rate v^cc/Va scaling as i5/L <C 1; this 
would not be rapid enough to be of any practical inter- 
est. Thus, we come to Conclusion I: irrespective of the 
actual microphysics inside the layer, Petschek's mecha- 
nism, or a variation thereof, is absolutely necessary for 
sufficiently fast astrophysical reconnection, including so- 
lar flares (and even micro- and nanoflares). 

Because Petschek's model was able to reproduce the 
very short observed flare timescales, it quickly became 
very popular and people have believed in it for the next 
20 years. With the advent of computer simulations, how- 
ever, the original Petschek model came under flrc. In 
particular, several numerical studies (e.g., Biskamp 1986; 
Scholer 1989; Ugai 1992, 1999; Ma & Bhattacharjee 
1996; Uzdensky & Kulsrud 1998, 2000; Erkaev et al. 
2000, 2001; Biskamp & Schwarz 2001; Birn et al. 2001; 
Malyshkin et al. 2005; Cassak 2005, 2006) showed that 
in resistive MHD with uniform resistivity (and, by infer- 
ence, with resistivity that is a smooth function of plasma 
parameters, e.g., Spitzer, see Biskamp & Schwarz 2001) 
Petschek's mechanism fails and Sweet-Parker scaling ap- 
plies instead. The basic physical reason for this, as has 
been elucidated analytically by Kulsrud (2001), (see also 
Uzdensky & Kulsrud 2000 and Malyshkin et al. 2005) 
is that the transverse (i.e., perpendicular to the current 
sheet) magnetic field component, which is necessary to 
sustain Petschek's shocks, is swept away by the Alfvenic 
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flow out of the layer so rapidly that it cannot be regen- 
erated by the nonuniform resistive merging of the recon- 
necting magnetic field component. 

In addition, a strong evidence for the existence 
of a slow Sweet-Parker reconnection mode in col- 
lisional resistive-MHD plasmas has been obtained 
from recent laboratory studies in the Magnetic Re- 
connection Experiment (MRX; Yamada et al. 1997) 
in the high-coUisionality regime (Ji et al. 1998; 
Trintchouk et al. 2003; Kuritsin et al. 2006). 

This enables us to draw Conclusion II: In the col- 
lisional regime, when classical resistive MHD applies, 
Petschek's fast reconnection mechanism does not work; 
slow Sweet-Parker reconnection process takes place in- 
stead. 

A very natural question to ask next is whether fast 
Petschek-like reconnection possible in a collisionless 
plasma, where resistive MHD does not apply. There is 
now a growing consensus that the answer to this question 
is YES. In Solar and Space physics, of course, there has 
long been plentiful observational evidence for fast colli- 
sionless reconnection, in solar flares (e.g., Tsuneta et al. 
1984, 1992, 1996; Masudaet al. 1994, 1995; Shibata 1995, 
1996; Yokoyama et al. 2001), and in the Earth magne- 
tosphere, both in the magnetopause (e.g., Mozer et al. 
2002, 2003, 2004) and in the magnetotail (Oieroset et al. 
2001; Nagai et al. 2001). The evidence for fast collision- 
less reconnection has been further significantly strength- 
ened by recent laboratory measurements in the MRX 
(Ji et al. 1998, 2004; Yamada et al. 2006) and other de- 
vices (e.g.. Brown 1999; Egedal et al. 2007; Brown et al. 
2006). These laboratory measurements, however, have 
not been able to elucidate the special role of the Petschek 
mechanism in accelerating reconnection. On the other 
hand, over the past decade or so, several theoretical 
and numerical studies have indicated that fast reconnec- 
tion enhanced by the Petschek mechanism (or a vari- 
ation thereof) does indeed take place in the collision- 
less regime. Moreover, it appears that there may be 
two mechanisms of collisionless reconnection. Physically, 
these two mechanisms are very different from each other; 
nevertheless, they both appear to lead to the establish- 
ment of a Petschek-like configuration, which enhances 
the reconnection rate. The two regimes in question are: 

• Hall-MHD reconnection, involving two-fluid effects 
in a laminar flow configuration (e.g., Mandt et al. 
1994; Kleva et al. 1995; Biskamp et al. 1995, 
1997; Ma & Bhattacharjee 1996; Lottermoser & 
Scholer 1997; Shay et al. 1998, 1999, 2001; Hesse 
et al. 1999; Bhattacharjee et al. 1999, 2001; 
Birn et al. 2001; Rogers et al. 2001; Pritchctt 2001; 
Breslau & Jardin 2003; Cassak et al. 2005, 2006; 
Daughton et al. 2006). 

• Spatially-localized anomalous resistivity due to 
plasma micro-turbulence that is excited when the 
current density exceeds a certain threshold (e.g., 
Coppi & Friedland 1971; Smith & Priest 1972; 
Coroniti & Eviatar 1977; Kulsrud 2001, 2005; Kul- 
srud et al. 2005). This seems to lead to a Petschek 
configuration with the length of the inner diffusion 
region being on the order of the resistivity localiza- 
tion scale (e.g., Ugai & Tsuda 1977; Sato & Hayashi 



1979; Ugai 1986, 1992, 1999; Scholer 1989; Biskamp 
& Schwarz 2001; Erkaev et al. 2001; Kulsrud 2001, 
2005; Uzdensky 2003; Malyshkin et al. 2005). The 
corresponding reconnection rate is then fast enough 
to explain the relatively-short timescales observed 
in solar flares (e.g., Kulsrud 2001, 2005; Uzden- 
sky 2003). 

Recent experimental evidence from the MRX exper- 
iment indicates that both anomalous resistivity (e.g, 
Ji et al. 2004) and HaU effect (Ren et al. 2005; Ya- 
mada et al. 2006; Brown et al. 2006; Egedal et al. 2007) 
are present and may be important in collisionless recon- 
nection. However, at present, it is still not clear which 
one of these two mechanisms works in a given physical 
situation (if at all). Also not known is whether these 
two mechanisms can operate simultaneously in a given 
system and how they interact with each other. 

In any case, from the point of view of our present 
discussion, the important thing is that one does get a 
Petschek-enhanced fast reconnection if the plasma is col- 
lisionless (in the sense that will be specified below, see 
§ 12. 2p . Thus, we can draw Conclusion III: a Petschek- 
enhanced fast reconnection does happen in the collision- 
less regime. 

To sum up, there are two regimes of astrophysical mag- 
netic reconnection: slow Sweet-Parker reconnection in 
resistive-MHD with classical collisional resistivity, and 
fast Petschek-like reconnection in collisionless plasmas. 
In other words, in order for reconnection to be fast, it 
needs Petschek's mechanism to operate and that in turn 
requires the reconnection layer to be collisionless. 

In fact, whenever we observe violent and rapid ener- 
getic phenomena that we interpret as reconnection, it is 
always in relatively tenuous plasmas. I am not aware of 
any counter-examples and would be would be very inter- 
ested in learning about them. Is there any evidence for 
fast large-scale reconnection events in collisional astro- 
physical environments? 

2.2. The Fast Reconnection Condition: Zero Guide 
Field Case 

How can one quantify the transition between the two 
regimes? For clarity, let us first consider the case with 
zero guide field, — 0. [In the ^ case some of 
the arguments and results presented in this section will 
be somewhat modified; however, the main concepts and 
conclusions will remain similar, as we show in the next 
section (see ? l2.4p .] 

When there is no guide field (strictly anti-parallel field 
merging), the condition for fast reconnection can be for- 
mulated roughly as (Kulsrud 2001, 2005; Uzdensky 2003; 
Cassak et al. 2005; Yamada et al. 2006) 

SsP <d, = —. (2) 

What this condition means is the following. As a recon- 
nection layer forms, its thickness S becomes smaller and 
smaller. If condition ^ is not satisfied, then this thin- 
ning saturates at ^ = Sgp', then, reconnection proceeds 
in the slow Sweet-Parker regime. On the other hand, if 
condition ^ is satisfied, then various two-fluid and/or 
kinetic effects kick in as soon as S drops down to about di 
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or so, i.e., well before the resistive term becomes impor- 
tant. Then, the reconnection processes necessarily in- 
volves coUisionless, non-classical-resistive-MHD physics. 
This results in a rapid increase in the reconnection rate, 
as has been recently documented experimentally in the 
MRX (Yamada et al. 2006). 

Note that, by construction, Sgp that enters in equa- 
tion ([2]) is to be calculated using the classical Spitzer 
resistivity. Therefore, and since we are talking about col- 
lisional and coUisionless reconnection, it is instructive to 
express the above condition in terms of the classical elec- 
tron mean free path due to Coulomb coUisions, Ag^mfp- It 
is pretty straight-forward to show (Yamada et al. 2006) 
that 

^ - V'7/3e ^ (3) 

V ^e,mfp J \ rrii J 

Here, f3e is the ratio of the electron pressure inside the 
layer to the pressure of the reconnecting magnetic field 
component (Bq/Stt) outside the layer. Thus, using equa- 
tion ([2]), we see that a reconnection layer is coUisionless 
when 

L < Lc = Ae,mfp y/m~JP~me ~ 40 /J^^/^ Ae^,„fp . (4) 

The coUisionless reconnection condition in a similar form 
has been first published by Yamada et al. (2006). Their 
condition actually differs from equation (j4]) by a factor 
of 2. Since the discussion here is very qualitative, I re- 
gard this difference as unessential. Moreover, to make 
the present discussion more clear, I systematically ignore 
numerical factors of order f everywhere in this paper. 

We can go a little bit further. The mean free path 
Ae,mfp can be written as 



A, 



e,mfp 



7- lO^cmn^fo^Ty^ 



(5) 



nge"* log A 

where we set the Coulomb logarithm equal to 20 and 
where nio and T-j are the central layer electron den- 
sity Ue and temperature Te given in units of 10^° cm~'^ 
and 10^ K, respectively. Combining equations ^ 
and ([5]), the criterion for fast coUisionless reconnection 
can now be formulated as a condition on the layer's 
length L in terms of the values of rie and Tg at the center 
of the Sweet-Parker reconnection layer: 

L < Le = 40 Ae,„>fp ~ 3 • lO^cm p-^/^ n'^ . (6) 

Notice that the critical layer length Lc strongly de- 
pends on the central electron temperature Tg. The 
main cause for this high sensitivity is the strong tem- 
perature dependence of the electron mean free path and 
hence Spitzer resistivity. It is therefore very important 
to figure out what the electron temperature inside the 
Sweet-Parker reconnection layer should be. In a situa- 
tion where the ambient (upstream) pressure is already 
high, /3upstroam > 1, the conversion of magnetic en- 
ergy to plasma via reconnection is not going to affect 
the plasma temperature and density considerably. In 
this case, both the plasma temperature and the density 
should be roughly uniform across the layer; then one can 
just substitute their upstream values into equation 
The story then ends here. 

However, in this paper we will be mostly inter- 
ested in environments that are (globally) magnetically- 
dominated (and hence almost force- free), such as the so- 
lar corona and coronae of other stars and accretion disks. 



In this case, the ambient thermal pressure is smaller than 
the reconnecting magnetic field pressure. It is then im- 
portant to distinguish the temperature inside the layer 
from the ambient coronal temperature (which is much 
lower) . The simplest first step towards determining is 
to note that, in the case with zero guide field, the condi- 
tion of pressure balance across the layer dictates that 
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^0 



1. 



(7) 



Next, assuming that Tg inside a Sweet-Parker reconnec- 
tion layer is approximately equal to (and in any case not 
much smaller than) the ion temperature Ti — the self- 
consistency of this assumption will be demonstrated at 
the end of this section — we see that f3e should generally 
be close to 1/2 (and in any case of order 1). 

Then, using the pressure-balance condition ([7]), we can 
express the central electron temperature in terms of Bq 
and Ue as 



Te = 



2kBne 



lA-10'KBl^n^o\ 



(8) 



where Bi ^ is the outside magnetic field Bq expressed 
in units of 30 G. Upon substituting this expression into 
equation ([5]), we get 



Ae,mfp ^ 1.5- lO^cmnfo'^L. (9) 

and so the fast coUisionless reconnection condition can 
be written in the final form as 

L < Lc{7i, So) ~ 6 • lO^cm nlg^ B^ ^ . (10) 

We see that the condition L < is often satisfied in the 
solar corona (e.g., in solar flares) and is definitely always 
satisfied in the Earth magnetosphere. 

2.3. Plasma Temperature and Density inside the 
Sweet-Parker Reconnection Layer 

Note that the density that enters in the above expres- 
sions is, by definition, the density at the center of the 
reconnection layer. In general, it is not known a priori 
and we would like, instead, to get expressions involving 
the background (upstream) density. Equation ([7]) gives 
us one relationship between the two desired quantities 
(central rig and Tg). However, we need to know both of 
them and by itself, equation ([7]) does not tell us whether 
the central pressure is increased (from a very low ambi- 
ent value) to the required equilibrium level by raising the 
density or the temperature. I will now argue, however, 
that, in the regime that is relevant here, the latter is 
more likely to be the case, i.e., that the pressure increase 
is mostly due to the increased temperature, whereas the 
density should change relatively little. 

Indeed, as is evident from the above discussion, the 
main reason for the strong density dependence in equa- 
tion pop is the Tg-dependence of the classical Spitzer 
resistivity, combined with the pressure balance condi- 
tion ([7]). Therefore, what we are mostly interested in here 
are the central electron density and temperature in the 
context of the resistive Sweet-Parker theory. To deter- 
mine them, we need more detailed information. Namely, 
we need to consider one more equation that we have not 
discussed up until now — the equation of energy conser- 
vation. In general this equation should take into account 
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Ohmic heating and all the relevant loss mechanisms, such 
as advection, radiation, and electron thermal conduction. 

Let us consider a fluid element as it travels through the 
reconnection layer. During its transit, the fluid element is 
subject to Ohmic heating and to all the above loss mech- 
anisms. The characteristic time for the advection term 
is of course just the Alfven transit time ta ^ L/Va — 
the time a fluid element spends inside the layer. During 
this time the element is continuously heated by Ohmic 
heating, whose rate per unit volume can be estimated 
within the Sweet-Parker theory as 

" 4tt (52 
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where the magnetic diffusivity rj is related to the resistiv- 
ity rj' via r] = rj'c^ /An.^ Thus, during its transit, the fluid 
element acquires just the right amount of heat to raise 
its temperature to about the equipartition level given by 
equation ([8]) (Uzdensky 2003). We now need to see under 
what conditions one can ignore various loss mechanisms, 
namely, radiation and thermal conduction. 

First, it is easy to see that, in the case of the solar 
corona, the radiative losses are indeed negligible on the 
time scales of interest (ta)- Indeed, the radiative cooling 
time is determined by the cooling function Q{T) via 

2cynkBT _ "iksT 
""■■"^ ^ n^Q{T) ~ ^^Q{T) ■ ^ ' 

In the temperature range of interest to the solar corona, 
10^ K < Te < 10^ K, the cooling function varies between 
10~23 and 10~22 gj-g sec~^ cm^ (e.g., Rosner et al. 1978; 
Priest 1984, p. 88; Cook et al. 1989). Then, taking 
n = 10^° cm~^, the typical radiative cooling times are of 
order 10 min (for T — 2 - 10® K) and longer (up to many 
days for T = 10^ K). In any case, this is much longer 
than the characteristic Alfven crossing time ta which is 
usually a few seconds. 

On the other hand, the smallness of the radiative losses 
in other astrophysical situations cannot be taken for 
granted, especially in High-Energy Astrophysics appli- 
cations, such as coronae of accretion disks around black 
holes and neutron stars and also gamma-ray bursts. If 
radiative losses are important, then the above estimates 
for the central temperature would no longer apply and 
correspondingly the condition for transition to the fast 
reconnection regime would have to be modified. 

Next, going back to the solar corona, we need to con- 
sider the effect of electron thermal conduction on the 
central temperature. Note that electrons are still well 
magnetized throughout most of the reconnection layer 
(this will be even more so in the presence of a guide 
field). Therefore, the electron thermal conduction is 
highly anisotropic, and so we shall estimate the effects 
of parallel and perpendicular thermal conduction sepa- 
rately. We shall start with the heat conduction along the 
magnetic field. Since we interested in the coUisional case, 
Ae.mfp ^ L, the parallel heat transport is in the diffu- 
sive regime and can be described by a one-dimensional 

^ I would like to add one extra word of caution here. The Sweet- 
Parker theory that we are relying on here does not take into account 
the large variation of the Spitzer resistivity across the layer as the 
electron temperature increases by a large factor with respect to its 
relatively low ambient value. This brings in an extra uncertainty 
into our estimates. 



random walk of electrons along the magnetic field. The 
characteristic parallel thermal conduction time-scale can 
then be estimated as 



^e.cond, || 



A, 



e,mfp 
Ve,th 



e,mfp 




where ta and f3 correspond to the reconnecting field Bq. 
Using expression we can write this expression as 



^e.cond. 



(14) 



In the magnetically-dominated environment of the so- 
lar corona, where the thermal pressure outside of the 
reconnection region is small, we expect /3 ~ 1; it can 
conceivably be smaller, but we never expect it to be sig- 
nificantly higher than 1. Therefore, we estimate that 



tld,|| > TA 



L 



(15) 



Thus, we see that if we are in the regime of inter- 
est — that is, in the expected regime of applicability 
of the coUisional Sweet-Parker theory, L > Lc, — then 
''e.cond.ll > TA and hcncc the energy losses due to parallel 
electron thermal conduction are not important. 

A similar line of arguments can be made to show that 
the characteristic time for perpendicular (across the re- 
connection layer) electron thermal conduction is in fact 
automatically always comparable to the Alfven cross- 
ing time. Indeed, the cross-field electron conduction 
timescale in the Sweet-Parker regime can be estimated 



as 



'^e,cond,_L 



D. 



e.± 



(16) 



where De,j_ ~ pl'^e.th/ ^e,mip is the classical electron 
coUisional diffusivity across the magnetic field. Using 
/?e ~ 1, we see that re_cond,_L should scale as 



'^e,cond,_L 




(17) 



Upon plugging equation ([3]) into this expression, we im- 
mediately find that all of the factors on the right-hand- 
side cancel and we simply see that 



T£,cond,± ^ Ta 



(18) 



i.e., that the two timescales are automatically compara- 
ble to each other. 

Combined with our previous result regarding the ef- 
fect of the parallel thermal conduction, we conclude 
that energy losses due to both parallel and perpendic- 
ular electron thermal conduction out of the layer are, at 
best, only marginally important in the coUisional Sweet- 
Parker regime. Heat losses are still present but are not 
likely to cause a decrease in the central electron temper- 
ature by more than a factor of order unity. This means 
that the jump in the plasma pressure, required to main- 
tain the pressure balance with the outside reconnecting 
magnetic field, should be attributed mostly to the in- 
crease in the plasma temperature (due to Ohmic heat- 
ing) , whereas the density should not vary across the layer 
by more than a factor of a few. Thus, for our rough es- 
timates, the density that enters in our estimate ^ for 
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the central temperature can be taken to be the ambi- 
ent plasma density. As long as we are in the collisional 
regime, L > Lc, this result is consistent with the energy 
balance inside the layer that includes Ohmic heating, 
heat advection by the bulk plasma outflow, and electron 
thermal transport. As we shall see in the next section, 
this will be true even in the presence of a guide field. 

It is also interesting to check that the assumption Ti ~ 
Tf, is not strongly violated. In principle, this could 
be a worry, since the collisional electron-ion energy- 
equilibration rate is suppressed by the large mass ratio; 
hence, in general, the electron and ion temperatures in 
the layer need not be equal. For example, ions might have 
been much hotter than the electrons and hence might 
have provided the bulk of the pressure support against 
the outside magnetic field. The electron temperature 
in this case would have been far below the equiparti- 
tion value (about 10^ K). Correspondingly, the electron 
mean-free path would be much lower than that given by 
equation Q. 

The electron- ion temperature equilibration time, teq, 
is by a factor of rrii/me longer than the electron colli- 
sion time, Tg. Making use of the the pressure balance 
Cg ~ (Te/m^)^/^ ^ Va, we can express in terms 
of Ae^infp and the Alfven speed as ~ Ae,mfp/''^e,th ~ 

(Ae,mfp/cs) ^m^jmi ~ (Ae,mfp/VA) \fmjmi. Then, we 
can compare teq with the Alfven crossing time, which 
is the characteristic time a a fluid element spends inside 
the layer: 



Teq ~ TA 



TA 



(19) 



L y me L 
Thus we see that, if we are in the collisional regime as 
defined by equation (jl]), the electrons and ions experi- 
ence enough collisions with each other to equalize their 
temperatures while they transit through the layer. Thus, 
Ti ~ Te should be a decent approximation in the Sweet- 
Parker regime. 

Note that the plasma density also enters the collision- 
less reconnection condition via di. The straight-forward 
logic of the above picture dictates that this quantity is 
to be estimated in the collisional regime. However, it 
is plausible that a corresponding estimate based on the 
collisionless regime also has some relevance; in this case 
it is interesting to compare the two. The task of estimat- 
ing the central rig in the collisionless regime is more diffi- 
cult than estimating it in the Sweet-Parker regime, since 
there are many uncertainties here. However, it is also less 
important, since most of the strong density dependence 
in equation (|10[) comes from the the steep temperature 
(and hence, indirectly, density) dependence of the Spitzer 
resistivity. Nevertheless, we still shall try to discuss the 
question of whether the density can vary significantly 
across a collisionless reconnection layer. If the fast colli- 
sionless reconnection is due to anomalous resistivity, then 
the effective Ae,mfp is determined by wave-particle colli- 
sions and hence will be relatively short. This will result, 
again, in suppression of parallel electron thermal conduc- 
tion and hence in high temperature and correspondingly, 
according to the pressure balance condition ([7]) , in a rela- 
tively mild density change with respect to the upstream 
value (see also Uzdensky 2003). Then the above esti- 
mates are likely to be still valid (although in this case 
one might have to deal with a possibility that Ti 3> Tg). 



On the other hand, if the fast reconnection is due to 
the Hall effect, then we also may expect the density to 
be roughly uniform across the layer. This is because 
the ions become demagnetized on scales less than di and 
hence cannot develop structures on smaller scales. The 
dynamics of ions inside the reconnection layer can be 
summarized by noting that they are pulled directly into 
the layer by the bipolar electric field, which accelerates 
them up to the Alfven speed, and then are ejected out 
along the layer by the Lorentz force associated with the 
quadrupole out-of-plane magnetic field (e.g., Uzdensky & 
Kulsrud 2006). At any rate, one does not expect to see 
a significant variation of ion, and hence electron, den- 
sity across the Hall reconnection layer (apart from the 
di-thick density depletion layers that run along the sep- 
aratrices, as reported by Shay et al. 2001). 

I also would like to make a little side remark about my 
use of terms here. In this paper I use the terms "colli- 
sional" and "collisionless" a lot and I would like to clarify 
what I mean by them. Notice that the collisionless re- 
connection condition ^ is different from the customary 
usage of the term "collisionless". Often, by "collision- 
less" one means a regime in which fie > v^^i, i.e., when 
the electron gyro- frequency is faster than the electron-ion 
collision rate. Equivalently, this means that pe < Ae,mfp- 
In this paper, on the other hand, the terms "collisional" 
and "collisionless" involve comparing Ae,mfp not with the 
electron gyro-radius, but instead with the system size L 
divided by the square root of the ion/electron mass ratio, 
as seen from 

Finally, I would like to make a comment about the use 
of the term "collisionless" . As I mentioned in § [U in 
the present paper this term is understood in the sense of 
equation ([4]). That is, a system is called "collisionless" 
when Ae.mfp ^ ^/40, roughly speaking. This definition 
is different from the condition Ag^mfp > Pe that is of- 
ten used. Note that, when Ae,mfp > Pe, or, equivalently, 
^^e > i^e.i, then the Hall term in the generalized Ohm law 
formally cannot be neglected. In plasmas of interest to 
this paper, e.g., in the solar corona, this condition is al- 
ways very well satisfied. In this sense, the coronal plasma 
is definitely always collisionless. Our collisionless recon- 
nection condition ([2|), on the other hand, is not always 
satisfied. In fact, there is a large and astrophysically- 
important region in parameter space where (5sp > di but 
where the plasma is nevertheless collisionless in the sense 
that > Vf.^i, and where therefore the Hall effect may be 
significant (see also Cassak et al. 2005). In the above dis- 
cussion I have assumed that, in this intermediate regime, 
reconnection still proceeds slowly, in the purely resistive 
Sweet-Parker mode, as suggested by numerous resistive- 
MHD simulations. One should be aware, however, that 
most of these simulations completely ignore the Hall ef- 
fect; thus, it is not clear to what extent their results can 
be relied upon in this intermediate regime. In fact, the 
recent numerical work by Cassak et al. (2005), who in- 
clude both resistivity and Hall effect, indicates that both 
slow Sweet-Parker and fast Hall modes may operate in 
the intermediate regime, depending on the previous his- 
tory of the system. In particular, they observe that, when 
one considers a gradual initial thinning of the reconnec- 
tion layer, the system always stays in the slow regime, 
even if fie > t'e.i, as long as (Ssp > di. Since it is ex- 
actly this scenario that is of interest in the context of 
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our model, we feel we can use condition ([2]) with confi- 
dence. 

2.4. Reconnection with a Guide Field 

In the real world the zero-guide-field case considered 
in the previous section is rare. It may be encountered in 
some special circumstances; for example it may be rele- 
vant to reconnection in the Earth magnetotail. Gener- 
ally, however, there is some guide field, 0, present. 
This is especially so in the context of the nanoflare model 
for heating of solar coronal loops by braiding of individ- 
ual elemental flux strands that comprise a larger loop. 
In this case the guide field (i.e., the axial field along the 
general loop direction) is always present and is generally 
much stronger than the reconnecting field Bq. Another 
situation where guide field is dominant over the recon- 
necting magnetic field is laboratory plasma devices, such 
as tokamaks. Note that reconnection releases the energy 
of the weaker i?o-field; the energy of the guide field is 
not available. However, the presence of a guide field, 
especially if it is strong, may, in principle, change the 
structure of the reconnection layer. Therefore, in this 
section we consider the effect of the guide field on our 
collisionless reconnection condition. 

The presence of the guide field modifies the collision- 
less reconnection condition. This is because, whereas the 
characteristic layer thickness in resistive MHD is not af- 
fected by the guide field and is still given by the same 
Sweet-Parker scaling ^ , the characteristic scale for col- 
lisionless reconnection layer changes. In particular, ac- 
cording to the published literature (e.g., Kleva et al. 
1995; Rogers et al. 2001; Cassak et al. 2007), the layer 
thickness for Hall collisionless reconnection becomes the 
ion acoustic Larmor radius, ps ^ Cs/^i, where fli corre- 
sponds to the total upstream magnetic field and where 
Cs ~ {Te/rriiY/'^ is the ion-acoustic speed. This result 
has been recently confirmed experimentally (Egedal et al. 
2007). Then, the criterion for the onset of fast reconnec- 
tion in the Hall-MHD regime becomes (e.g., Cassak et al. 
2007) 

^SP < Ps . (20) 

According to the general scheme presented in the pre- 
ceding section, we first need to consider the Sweet-Parker 
layer. Again, we want to estimate the value of the Spitzer 
resistivity at the center of the reconnection layer. There- 
fore, we need to determine the central electron temper- 
ature. One of the difficulties here is that, if there is a 
guide field, then the cross-layer force-balance in the form 
of equation ([7]) is no longer applicable. This is because 
the pressure of the reconnecting component of the out- 
side magnetic field, Bq/8tt, no longer has to be balanced 
at the center of the layer by the plasma pressure alone. It 
can be partially balanced by an increase in the guide field 
pressure due to a modest compression. In particular, a 
very strong guide field, Bz S> -Bq: dominates the pressure 
both inside and outside of the current layer and can be 
changed only slightly. Therefore, it essentially ensures 



incompressibility : 



The temperature then 



decouples from the pressure balance and has to be deter- 
mined from some other considerations, namely, from the 
energy equation. Fortunately, most of the analysis pre- 
sented in the previous section can be carried over to the 
strong-guide-field case without much change. In partic- 



ular, the plasma experiences Joule heating as it travels 
through the layer and, as long as one is in the collisional 
regime, the parallel electron thermal conduction is not 
able to cool the layer by a large factor. Therefore, we ex- 
pect the temperature inside the layer to increase roughly 
to the equipartition level given by equation ([8]). However, 
this equation can no longer be regarded as equality, but 
just as a rough estimate. For the purposes of this paper, 
though, this is good enough. 

Now, let us discuss the right-hand side of condi- 
tion ((20|) . If the guide field is not very large compared 
with the reconnecting field, then one expects ps ^ di, 
and hence the criterion ((2|) is not changed substantially. 
However, in the strong-guide-ficld case, we expect to have 
Ps di (assuming, as we always do throughout this pa- 
per, that upstream /3 <^ 1). Then, the criterion for tran- 
sition to fast collisionless Hall reconnection is modified 
structurally. It now becomes 



<5sp < Ps ^ di 



1/2 



d. f , (21) 



where /?e is based on the central electron temperature 
and density and on the upstream reconnecting field com- 
ponent Bq. Note that, according to logic of our paper, 
this /3e is to be estimated in the collisional regime. Nev- 
ertheless, it may also be useful to get some feeling for 
what /3e (and hence p^) should be immediately after the 
transition to fast reconnection. As we noted above, a 
strong guide field makes the plasma nearly incompress- 
ible, so that rig is roughly uniform. The question then 
is what one should take for the electron temperature in- 
side the layer? There is a considerable uncertainty here. 
In particular, there is no good reason for Ti to be close 
to Te in the collisionless regime. For example, Hsu et al. 
(2001) present experimental evidence for strong ion heat- 
ing Ti ^ Te inside a collisionless reconnection layer in the 
no-guide-field case in the MRX, although the rise in Ti 
is considerably smaller in the presence of a guide field. 
In any case, with all these reservations, the simplest ap- 
proach one can adopt is to take /3e 1 in the above 
expression. Then, the condition for fast collisionless re- 
connection in the strong guide field case becomes 



^SP 
Ps 



L 



1/2 / s 1/4 n 

' m,A B, 



nil 



Bo 



< 1, 



(22) 



which, after our usual manipulations, can be rewritten 
as 

2 



L < Lc{Bz) — \l Ae,infp ( 77^ ) 

V me \BzJ 



i6- 10^ cmn7^ B'^ 



Bn 



10 ^1.5 [^J (23) 

Thus, the critical global length Lc in this case becomes 
smaller by a factor [Bz/BoY and so the collisionless re- 
connection condition is harder to satisfy. Also note that, 
for fixed n,, and Bz, Lc becomes very sensitive to the 
reconnecting field component: Lc 



£.Q. 



2.5. Caveats and Alternatives 

With all the above said, there are still many caveats 
and alternative ideas that may potentially alter the above 
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picture. It is important to be aware of them and try to 
gage their effect on the existence and the form of the fast 
colhsionless reconnection condition. I wiU discuss some 
of them in this section. 

First, there are a number of ideas that may, if cor- 
rect, represent a potential chaUenge to our Conclusion II 
(that collisional resistive-MHD reconnection is slow) . For 
example, Lazarian & Vishniac (1999) suggested that re- 
connection may be relatively rapid in pure resistive MHD 
in the presence of externally imposed three-dimensional 
MHD turbulence (see, however, Kim & Diamond 2001 
for a dissenting view). As far as I know, this idea has 
not yet been tested in three-dimensional (3D) numeri- 
cal simulations, although it is definitely an interesting 
possibility that deserves further study. Such studies are 
indeed already underway (Vishniac 2007, private com- 
munication). Furthermore, even in the two-dimensional 
case, resistive-MHD Sweet-Parker-like current layer may 
become unstable to tearing instability when its aspect 
ratio exceeds a certain large number (probably a hun- 
dred) (Bulanov et al. 1978; Forbes & Priest 1983; 
Biskamp 1986; Lee & Fu 1986; Malyshkin 2005, private 
communication; Loureiro et al. 2007). Tearing may effec- 
tively result in an enhanced hyper-resistivity that would 
dominate over the normal resistivity (e.g., Strauss 1988). 
What will be the effect of the nonlinear development of 
this instability on the overall reconnection rate is still 
not clear. In addition, Dahlburg et al. (1992) identified a 
mechanism, based on a secondary ideal-MHD instability, 
that leads to the excitation of fully three-dimensional tur- 
bulence within a reconnecting current layer. Finally, in a 
recent numerical simulation with both uniform resistivity 
and Hall effect, Cassak et al. (2005) observed an interest- 
ing hysteresis phenomenon: as the resistivity was turned 
down, there was a transition from slow Sweet-Parker to 
fast Hall reconnection, in accordance with condition 0. 
However, as the resistivity was turned back up to a level 
higher than that corresponding to the system still 
stayed in the fast mode. This puzzling behavior needs to 
be investigated further. 

Second, there are recent studies relevant to Conclu- 
sion III (that colhsionless reconnection is fast). In par- 
ticular, even though most of the authors agree that Hall 
reconnection is much faster than Sweet-Parker, there is 
still a significant disagreement about exactly how fast it 
actually is. For example, Drake and his collaborators 
report a universal reconnection rate of about 0.1 Va, 
independent of the electron microphysics that is ul- 
timately responsible for breaking the field lines (e.g., 
Biskamp et al. 1997; Shay et al. 1998, 1999, 2001; Cas- 
sak et al. 2005; see also Birn et al. 2001). On the other 
hand, others argue that it should not be universal but 
should instead be a function of the system parameters, 
such as the electron/ion mass ratio (e.g., Bhattacharjee 
et al. 2001). In addition, Daughton et al. (2006) re- 
cently raised a concern about the validity of periodic 
boundary conditions used in the majority of previous 
numerical studies. In particular, they argued that the 
reported very high reconnection rate may be an arti- 
fact of such boundary conditions. Furthermore, they 
performed 2D particle simulations with the appropriate 
open boundary conditions and observed that, after an 
initial transient peak consistent with the high value re- 
ported by Shay et al. (1999, 2001), the reconnection 



rate went down and settled at a significantly lower level 
(Daughton et al. 2006). Similar conclusions were reached 
independently by Fujimoto (2006) using periodic 2D par- 
ticle simulations but in a very large. After a year-long 
debate, however, this controversy finally seems to have 
been resolved. It is now agreed by both sides that the 
electron dissipation region develops a two-scale structure 
with an extended high- velocity electron jet in the outflow 
direction; however, the dimensionless reconnection rate 
still remains high (uroc/VA ^ di/L), although perhaps 
not quite as high as previously believed (e.g.. Shay et al. 
2007; Karimabadi et al. 2007). 

Finally, I would like to mention that the picture pre- 
sented in the above two sections, implies certain physi- 
cal conditions and thus is not expected to be universally 
valid in all astrophysical environments. For example, in 
the case of weakly-ionized plasma in molecular clouds 
(e.g., in the context of star formation), reconnection dy- 
namics is strongly modified by recombination processes 
and may become significantly faster than Sweet-Parker 
even in the collisional regime (Zweibel 1989; Heitsch 
& Zweibel 2003). Also, in the case of reconnection 
of super-strong magnetic fields in magnetar magneto- 
spheres {Bq > 10^^ G), the density of the released mag- 
netic energy is so large that the pressure inside the layer 
is dominated by radiation. The resulting temperature 
is so high that a large number of electron-positron pairs 
is produced; this means that the number of particles is 
not conserved, and so the above estimates would have 
to be modified (see, e.g., Thompson 1994; Lyutikov & 
Uzdensky 2003; Uzdensky & MacFadyen 2006). 

Other issues that are worth-mentioning and that may 
substantially complicate the simple physical picture pre- 
sented in this paper are: impulsive, bursty reconnection 
(e.g., Bhattacharjee et al. 1999; Bhattacharjee 2001; 
Daughton et al. 2006, 2007, in preparation); and recon- 
nection in three dimensions (e.g., Longcope 1996). 

3. SOLAR CORONAL HEATING AS A SELF-REGULATED 
PROCESS 

3.1. General Picture 

The physical picture presented in this paper has a num- 
ber of interesting implications for the solar corona that 
I will discuss in this section. 

At the present, the theory of magnetic reconnection 
has not yet reached the level of maturity that would al- 
low us to make accurate predictions regarding when the 
transition between the slow and fast reconnection modes 
will occur. Moreover, there is still no consensus on the 
exact nature of the fast colhsionless reconnection mecha- 
nism and on the actual scalings of the reconnection rate 
in either fast or slow reconnection (see section § 12.51 for 
more discussion). However, what is important for us in 
this paper is just the sole fact of existence of such a tran- 
sition and that it has to do with the plasma collisionality. 
The transition doesn't actually have to be razor-sharp, 
and in reality there may be not two but more different 
regimes with different reconnection rates and hence more 
than one such transition. Also, the transition condition 
may depend on other physical parameters, such as the 
presence of a guide field (see § 12. 4p or, perhaps, the ef- 
fectiveness of radiative cooling, which we have omitted 
in this paper. Whatever the case, the existence of two 
reconnection regimes with a coUisionality-related transi- 



Solar Coronal Heating 



9 



tion between them is a very important fact and we would 
like to exploit it in various astrophysical applications. In 
this section, I will apply this condition [namely, eqs. ^ 
and (fTO|) ] to the problem of solar coronal heating. 

I propose that coronal heating is a self-regulating pro- 
cess that works to keep the corona marginally collision- 
less, in the sense of equations ^ and (|3|) (see also Uz- 
densky 2006, 2007). 

According to Parker's nanoflare theory (Parker 1972, 
1983, 1988), as long as the twisting and braiding of coro- 
nal loops by photospheric footpoint motions and flux 
emergence episodes keep producing thin current sheets 
in the corona (e.g., van Ballegooijen 1986), magnetic dis- 
sipation in these current sheets leads to continuous coro- 
nal heating. In reality, this heating, of course, is not 
uniformly distributed; it is instead strongly intermittent, 
localized in both space and time (Rosner et al. 1978; 
hereafter RTV78). This basic picture of coronal heating 
has been shown to work in several 3D MHD numerical 
simulations (e.g., Galsgaard & Nordlund 1996; Hendrix 
et al. 1996; Gudiksen & Nordlund 2002, 2005). 

Now, the overall, integrated heating density, i.e., the 
rate of magnetic dissipation per unit volume, depends on 
the reconnection rate in these sheets, e.g., on whether 
reconnection is fast (e.g., Petschek) or slow (Sweet- 
Parker). This is actually a rather subtle point. Indeed, 
in a steady state, the energy dissipated in the corona per 
unit time should be equal to the power pumped into the 
corona magnetically from the solar surface. And it is 
not obvious why and how the energy-pumping rate de- 
pends on what happens in the corona. For example, if 
the corona were enclosed in a fixed volume, then the en- 
ergy dissipation per unit volume would be fixed (as in 
driven MHD turbulence in a box). However, it is impor- 
tant to recognize that the coronal volume is not fixed! If, 
for example, reconnection were to suddenly slow down, 
then more energy would be pumped in than the corona 
could dissipate; then, in order to accommodate this ad- 
ditional free magnetic energy, the corona would respond 
by increasing its scale-height. That is, coronal magnetic 
structures would grow in height until, finally, the total 
dissipation in the corona became equal to the total input 
from the photosphere. Because of the increased volume, 
the magnetic dissipation per unit volume is decreased. 

In addition, as coronal magnetic structures grow in 
height, especially when H > Rq, the amount of energy 
pumped into the corona by the footpoint motions may go 
down. This is because the work done by footpoint mo- 
tions is proportional to Vfp-Bhor^vcit/47r. As the coronal 
structures grow in height without increasing their lateral 
size, the horizontal field, i?hor, decreases, whereas the 
vertical field component, -Bvoit, does not change. Corre- 
spondingly, the overall power pumped from the photo- 
sphere into the corona goes down. 

As follows from equation (|10|). the reconnection mode 
is determined by the global scale L of the reconnection 
layer and by the basic physical parameters characteriz- 
ing the plasma in the layer (i.e, rig, Tg, and Bq). The 
typical values of L and i?o are determined by the scale 
and strength of the magnetic structures emerging from 
the Sun and by the scale of footpoint motions (e.g., the 
meso-granular scale). Therefore, for the purposes of the 
present discussion, I will regard L and Bq as fixed and 
then ask what determines the coronal density and tem- 



perature. 

Following this line of reasoning, let us invert equa- 
tion (jlOp and view it as the condition for the plasma 
density. That is, let us introduce a scale-dependent crit- 
ical density, below which the reconnection process 
transitions from the slow coUisional Sweet-Parker regime 
to a fast collisionless regime: 

nc{B, <Bo)^2- IQi" cm'^ ^-Vs ^ (24) 

where Lg is the global reconnection layer length ex- 
pressed in units of lO'' km. In the strong-guide-field case, 
Bz ^ Bq, the corresponding expression becomes: 

n.{Bz » Bo) -2- 10^° cm-3 ^-1/3 |^|o ^ 

(25) 

It is interesting to note that these values are close to 
those observed in active solar corona. I suggest that this 
is not a pure coincidence. 

The reason for this is that there is an important feed- 
back between coronal energy release and the coronal den- 
sity. This feedback is due to the fact that the gas high 
in the corona actually comes from evaporation from the 
surface along the field lines that just underwent recon- 
nection. In fact, it is essential to consider the coronal 
heating not only as the process of increasing the coro- 
nal plasma's temperature, but also as that of increas- 
ing its density (e.g., Klimchuk 2006; Aschwanden et al. 
2007). Indeed, a part of the energy released in a re- 
connection event is rapidly conducted along the field by 
fast (perhaps, nonthermal) electrons to the surface and 
is deposited in a denser photospheric and chromospheric 
plasma. This leads to a localized heating at the foot- 
points of the post-reconnected magnetic loops and to 
a subsequent chromospheric evaporation along them, a 
well-documented phenomenon in solar corona (e.g., As- 
chwanden et al. 2001a, b and Klimchuk 2006; see, how- 
ever. Brown et al. 2000 and Aschwanden et al. 2007 
for an evidence that chromospheric evaporation occurs 
in response to chromospheric, rather than coronal, heat- 
ing events). As a result, these loops become filled with a 
dense and hot plasma. The density rises and may now ex- 
ceed ric- This will shut off any further reconnection (and 
hence heating) involving these post-reconnective loops 
until they again cool down, which occurs on a longer, 
radiative timescale. 

Let us consider a couple of examples of how this might 
work. First, let us consider the case when the plasma 
density is relatively high and so the radiative cooling 
time. Trad, is relatively short compared with the timescale 
on which the coronal magnetic structures change due to 
footpoint motions, Tfp; however, it is of course still long 
compared with the fast reconnection timescale. The op- 
posite limiting case will be discussed later in this section. 
Let us suppose that due to field-line twisting, a recon- 
necting structure is set up in the corona at i = 0. Let this 
structure be characterized by the current sheet length L 
and the reconnecting field component Bq. Furthermore, 
let us further suppose that, initially, the ambient density 
of the background plasma is higher than nc{L, Bq). Then 
the reconnection layer is coUisional and reconnection pro- 
ceeds very slowly, in the Sweet-Parker mode. That is, 
there is almost no reconnection at all and hence coronal 
heating is weak. The surrounding plasma gradually cools 



10 



D. A. Uzdensky 



by radiation (and also, in general, by thermal conduc- 
tion) and the pressure scale height gradually goes down. 
The gas gradually precipitates to the surface. Then the 
density of the plasma entering the layer decreases and at 
some point becomes lower than the critical density. The 
reconnection process then suddenly switches to the col- 
lisionless regime. Petschek-like fast reconnection ensues, 
and the rate of magnetic energy dissipation greatly in- 
creases. A flare commences. Some fraction of the energy 
released by reconnection is transported by the electron 
conduction along the reconnected field lines down to the 
base, where it is deposited in the much denser photo- 
spheric or chromospheric plasma. This in turn leads to 
a massive evaporation along the same field lines. As a 
result, the newly-reconnected loops are now populated 
with relatively dense and hot plasma. They cool down 
only slowly via conduction and radiation losses, keeping 
their relatively high density for an appreciable length of 
time. If, during this time, these loops become twisted or 
somehow get in contact with other loops, they are now 
not likely to reconnect rapidly, since their plasma density 
is above critical. This inhibits further coronal heating in 
the given region. In fact, we can speculate that for any 
further outbursts of coronal activity in the given region 
to occur, one has to wait for the gas in post-reconnective 
loops to cool down significantly, which occurs on a longer, 
radiative timescale. 

Thus we see that, although highly intermittent and in- 
homogeneous, the corona is working to keep itself roughly 
at about the height-dependent critical density given by 
equation ([M]) or (P5|) . Correspondingly, the background 
coronal temperature should be such that results in a 
density scale-height that is just large enough to popu- 
late the corona up to the critical density level at a given 
height. In this sense, coronal heating regulates itself (Uz- 
densky 2006, 2007). 

An important remark is that loop brightness in soft 
x-rays and UV does not necessarily mean that this loop 
is actively undergoing fast reconnection at the given mo- 
ment of time. Instead, it just reflects the fact that the 
loop has a high plasma density (since radiation intensity 
is proportional to n^). This implies an episode of strong 
evaporation from the surface in the recent past, presum- 
ably caused by a preceding fast-reconnection coronal en- 
ergy release in this loop. The nearby coronal regions 
that are darker have evidently lower density, and they 
would be ripe for reconnection from the point of view of 
the collisionless reconnection condition. However, they 
probably are filled with a magnetic field that does not 
have complex topology with current sheets, since the un- 
derlying reason for coronal activity is always the genera- 
tion of small coronal current sheets as a result of field- line 
braiding driven by a complicated pattern of the footpoint 
motions on the solar surface. 

The above radiative-cooling-dominated scenario was 
presented here first only because it brings out more 
clearly the role that plasma density plays in controlling 
the reconnection regime. It also highlights the role of 
chromospheric evaporation in controlling the plasma den- 
sity. However, there is an alternative version of this pic- 
ture corresponding to the opposite situation. Tip < Tj-ad- 
In this case one can repeat the same arguments as those 
given above but instead of the slow initial evolution of 
the plasma density due to gradual radiative cooling, one 



can invoke the slow evolution of the reconnecting mag- 
netic field strength caused by the motion of the loop foot- 
points (as was described by Cassak et al. 2006). In other 
words, one can regard the ambient plasma density n as 
fixed on the timescale of interest, whereas the reconnect- 
ing field component Bq gradually increases (for simplicity 
we also keep L fixed). Then, the current sheet gradually 
builds up, while staying in the Sweet-Parker mode, and 
then rapidly switches over to the fast collisionless mode 
once Bq has exceeded a certain threshold. This scenario 
is in fact closer in spirit to the original Parker nanofiare 
model (Parker 1983, 1988), and also to the discussion by 
Cassak et al. (2005, 2006). Instead of defining nc{L, Bq), 
one can just as well rewrite the fast reconnection con- 
dition in terms of a critical current-sheet strength, i.e., 
the critical reconnecting magnetic field component. Be, 
expressed as a function of L and n: 

Bc{n, L) 2^ 20 G lI^^ n%'^. (26) 

In the case of a strong guide field, B^ ^ Bq, the critical 
reconnecting field is 

B,{n, L, B,) ^ 30 G lI'^ n}/' bI'^ . (27) 

where B^^2 = BJ(IOOG). 

As I will show below, in the long run one can ex- 
pect some sort of statistical steady state with rfp ~ Tj-ad ■ 
Then, this second scenario is just as likely to take place 
as the first one. In fact, in reality they probably both 
occur alongside each other. 

I would like to emphasize that the phenomenon de- 
scribed in the first part of this paper (i.e., the transition 
from slow to fast collisionless reconnection) is not the 
whole coronal heating story, although, in my view, it is 
an important part of it. Conceptually, it is probably 
about one quarter of the whole story. The other impor- 
tant pieces of the puzzle are: (i) complex dynamic inter- 
action of a large number of coronal magnetic loops, ulti- 
mately responsible for the generation of current sheets 
in the corona by footpoint motions; this requires the 
development of an appropriate statistical description of 
loops (e.g., Uzdensky & Goodman 2007, in preparation); 
(ii) chromospheric evaporation following a reconnection 
event and motions of plasma along the loop; and (iii) ra- 
diative (and thermal-conduction) energy losses and the 
resulting precipitation of the material back to the so- 
lar surface. All these essential physics ingredients (or 
at least some effective representation of them) will even- 
tually have to be incorporated into a single comprehen- 
sive theoretical (probably numerical) model of the solar 
corona. 

It is important to note that each of these processes 
comes with its own characteristic timescale. The recon- 
nection time in the fast regime. Tree, is very short. For the 
purposes of our discussion, we shall regard it as instanta- 
neous. The characteristic time for material evaporation 
following a nanofiare is also relatively short. The two 
longer time-scales in the problem are the characteristic 
timescale for generation and growth of new current sheets 
by photospheric footpoint motions, rfp, and the radiative 
cooling timescale Trad- As we shall show below, the sys- 
tem tends to adjust itself to a statistical equilibrium in 
which these two timescales are comparable (of order one 
hour in the solar corona). Finally, the longest time scale 
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of all is the Sweet-Parker reconnection time. It is much 
longer than one hour, and is therefore irrelevant in our 
problem. That is, the energy release during the Sweet- 
Parker phase is not important. However, this does not 
mean that the existence of the slow Sweet-Parker regime 
is not important; it is in fact extremely important, since 
it allows for a strong current sheet to build up, accumu- 
lating a large amount of free magnetic energy. 

This discussion brings us to one of the most interesting 
issues in solar coronal heating — intermittency. The key 
questions here are: what is the distribution of (spatial, 
temporal, and energy) scales on which magnetic dissipa- 
tion occurs? what scales are responsible for most of the 
energy release? In my opinion, the coUisionless reconnec- 
tion condition is, at the end of the day, very important 
for answering these questions. In particular, I think this 
condition plays a key role in establishing the spectrum 
of the energy release events, e.g., the observed power- law 
distribution of flare energies. 

3.2. Long-term Evolution of a Coronal Loop 

To illustrate the above points, let us consider the fol- 
lowing model problem. One can raise the question of 
whether the self-regulating behavior of the nonlinear dy- 
namical system described above can exhibit an oscilla- 
tory behavior (Parker 2006, private communication), es- 
pecially in light of a hysteresis-like phenomenon reported 
by Cassak et al. (2005). Of course, in general this issue 
is difficult to assess within the present model, since it 
first needs to be incorporated in the overall multi-loop 
dynamics (e.g., Uzdensky & Goodman 2007, in prepara- 
tion); that is, the energy release in a given fast recon- 
nection event does not, strictly speaking, directly affect 
the reconnection process that produced it. However, it 
does lead both to relaxation of the magnetic stress and 
also, importantly, to the increase in density in the post- 
reconnected flux tube (or more precisely, in the elemen- 
tary magnetic thread within the larger flux tube). This 
will in turn temporarily inhibit any reconnection involv- 
ing this given field thread even if a new current sheet 
is formed. This inhibition will last for some time un- 
til either the plasma density goes down because of the 
radiative cooling or the reconnecting field component in 
the new current sheet grows to a larger strength. 

Let us consider a simple illustrative example that 
demonstrates this oscillatory behavior and the secular 
evolution of the system on a longer timescale. Consider 
a composite coronal fiux tube, initially filled with a rela- 
tively low plasma density n. For simplicity, let us imagine 
that that this tube consists of two separate elementary 
flux threads that are being wrapped around each other 
steadily by the sheared photospheric motions. Conse- 
quently, a current sheet gradually builds up inside the 
tube, in the spirit of the nanoflare model. Again for 
simplicity, let us assume that the length L of this cur- 
rent sheet and the axial magnetic field strength in 
the tube stay constant at all times. Then, as the ele- 
mental flux threads are wrapped around each other, the 
reconnecting field component, Bq, gradually increases. 
Because the plasma density is low, the radiative cooling 
time is relatively long and therefore we can regard the 
plasma density constant during this phase. On the other 
hand, also because the density is low, the critical field 
strength is also low, according to equation (|27p. That 



is, the current sheet will become coUisionless relatively 
early, at some small value of Bq. Correspondingly, the 
amount of free magnetic energy attributed to the current 
sheet at that moment will also be small and the resulting 
reconnection event will be relatively weak. After the re- 
connection event, the magnetic field becomes unwrapped 
again (almost completely, since the hysteresis behavior 
reported by Cassak et al. 2005, 2007 implies that fast 
reconnection proceeds to the end); then the footpoint- 
driven process of current-sheet build-up starts all over 
again. Now, if the plasma density in the tube under 
consideration were fixed at a constant small value Ue, 
then (since L and B^ are also kept fixed) all the sub- 
sequent fast reconnection events would be triggered at 
the same critical value of Bq = Bc{n, L, B^)- Corre- 
spondingly, they would all be characterized by a rel- 
atively small amount of energy released, on the order 
of VB'^^Ue, L, Bz)/8Tr, where V is the volume of the re- 
gion involved in the given reconnection event. At the 
same time, since weaker current sheets are easier to pro- 
duce, there would be following each other at a relatively 
rapid rate. 

However, because of the mass exchange between the 
surface of the Sun and the corona, the ambient density 
is not fixed a priori. In particular, if one starts with 
a small initial density, the first few reconnection events 
will be energetically small, but, nevertheless, they will 
eventually lead to an increase in the plasma density via 
chromospheric evaporation. This will make the critical 
magnetic field stronger and hence the amount of free 
energy stored between reconnection events and released 
through them larger. At the same time, such events will 
be more rare, since it takes longer to build up a stronger 
current sheet. Therefore, an increase in plasma density 
shifts the dominant energy-release scale towards larger 
energies. Thus, in this picture, the plasma density in 
the tube increased over time in a step-like manner: it 
increases by a certain amount in the aftermath of each 
reconnection event and stays roughly flat between them. 
The amount of plasma pumped into the corona as a re- 
sult of each event is roughly speaking proportional to the 
magnetic energy released in this event. The latter is in 
turn a monotonically-increasing function of density, ac- 
cording to equation ([27|) . That is, we have, keeping L 
and Bz fixed: 5n ~ B'^{n) ~ n. Thus, the relative in- 
crease is then independent of density; it can be roughly 
estimated as follows. Consider a loop of constant cross- 
sectional area L^ and length L|| > L; the loop's volume 
is then Vioop ^ i|| L^- Assuming that the reconnection 
layer extends all the way along the loop, the volume V 
of the region involved in the reconnection event is sim- 
ilar: V Vioop ~ L^- Then, the magnetic energy 
released in the event is of order i^roc ^ {Bq/Stt) L|| L^. 
Next, let us say that a certain fraction e < 1 of this 
energy goes towards filling the entire loop with plasma 
evaporated from the surface. The coefficient e depends 
on the portion of the released energy that goes to the 
energetic electrons and is transported by them to the 
photosphere, times the fraction of this deposited energy 
that goes to material evaporation (as opposed to that im- 
mediately radiated away from the footpoints). Overall, e 
describes the fraction of the fiare energy that eventually 
will be radiated by post-coronal loops as soft X-rays, UV, 
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etc., but not hard X-rays. Immediately after the evapo- 
ration episode, the energy eSrcc is divided between the 
increase in the gravitational energy, (5i?grav, and in the 
thermal energy, SEth — 3 Sn fcsTioopMoop, of the coronal- 
loop plasma (for simplicity of discussion, I dropped the 
contribution to dEth associated with the increase in the 
gas temperature; in reality, of course, both n and Tioop 
should increase together). If the height of the loop is 
larger than the density scale-height, then these two con- 
tributions are likely to be comparable, according to the 
virial theorem. However, in the case of a compact loop, 
whose height is less than the thermal scale-height, the 
thermal energy increase will be larger than In 
either case, we expect 6Eth ~ e^'rec, from which we im- 
mediately find 



1 



2e 



(28) 



where /3ioop = 167rnfcBT'ioop/i?c the characteristic 
plasma-beta based on the critical reconnecting field. 
Substituting the expression (|27)) for we get 

Aoop — 0.07 nio Tioop. 6 i 5 0.07 Tioop.e Lg , 

(29) 

(30) 



and so 



in<:T-i r2/3 rl/3 



There are also numerical factors of order unity which we 
have ignored. It is reasonable to expect that this ratio 
should be roughly of order unity. For now let us just 
assume that it is not large, i.e., Sn/n < 1. 

As mentioned above, in a realistic scenario the loop 
temperature should also increase at each step, in concert 
with the density. This means that, to get into the coronal 
loop, chromospheric gas would need to be heated to a 
temperature that is higher than at the previous step. 
Then, each erg of the released energy would be able to lift 
a smaller amount of gas, i.e., the evaporation efficiency of 
coronal heating would decrease over time. To illustrate 
this point, let us, for example, adopt the famous RTV 
scaling for the loop temperature (RTV78), even though it 
was derived assuming steady heating and thus is not valid 
for the impulsive situation considered here (see below 
for more discussion). Recasting the RTV78 scaling in 
terms of the electron density instead of the pressure, we 
write T^^^ ~ 3 • 10^ K nj^^ Substituting this into 

equation (PO)) . we get 



dn 

— ^ o e n 



-1/2^2/3 j-l/3 j-\l2 
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'1/2 



(31) 



Next, the time separation between the steps also in- 
creases with the increased threshold magnetic field. In 
particular, let us assume that the footpoint driving twists 
up the field lines and thus creates and then strength- 
ens current sheets (or perhaps, quasi-separatrix lay- 
ers, see, e.g., Galsgaard & Nordlund 1996; Titov & 
Hornig 2002) at some constant rate. This is, of course, 
an over-simplification. In reality, current sheets may 
be created as a consequence of the nonlinear develop- 
ment of kink-like instabilities of twisted magnetic loops 
(Dahlburg et al. 2005). Then, the growth of Bo with 
time is not likely to be a smooth function, but instead 
may involve sudden jumps. Here, however, we would like 



to avoid such complications. Thus, let us consider, as an 
illustration, a simple model in which the reconnecting 
field component increases linearly between reconnection 
events, 

^^.a, (32) 

where 7 = const. Then, using equation (|27p . the time bt 
it takes the reconnecting field to grow a given critical 
value Bc{n) (after a complete relaxation during the pre- 
ceding reconnection event) is 5t = ^^^Bc{n)/Bz ^ n^/^. 
Therefore, as long as the relative increase in density at 
each step is not large, the long-term {t ^ St) evolution 
of the system can approximately be described by a dif- 
ferential equation: 



dn 
'dt 



Sn{n) 
5t{n) 



const , 



(33) 



corresponding to a linear rise n{t) ~ t. Accordingly, the 
emission measure of the loop should increase as . [Al- 
ternatively, if we neglect the evolution of the loop tem- 
perature in the above analysis and instead keep Tioop — 
constant, we then get 5n ^ n and hence nit) ~ t^, with 
the emission measure rising as i^.] 

This growth will continue until one of the following two 
processes will intervene. First, it may happen, as the 
density builds up, the critical value of Bq will become so 
large (e.g., a sizable fraction of S^), that the equilibrium 
shape of the entire loop will be affected. For example, 
the loop may become external-kink-unstable and become 
helical (sigmoidal loop), which with further twisting may 
result in a large-scale disruption with a catastrophic en- 
ergy release, a large flare. This scenario, by the way, sug- 
gests that coronal loops should gradually (on the time 
scale of hours) brighten up before a major disruption, 
because a higher plasma density (manifested as a higher 
emission measure) leads to a higher critical field strength 
and hence allows a larger amount of free magnetic energy 
to be accumulated without being prematurely dissipated 
via smaller reconnection events. 

On the other hand, there is another, less violent out- 
come that is also possible. Indeed, so far we have been ne- 
glecting radiative cooling. However, as the plasma den- 
sity gradually builds up due to a series of chromospheric 
evaporations caused by reconnection events, the radia- 
tive cooling time decreases (as n~^). At the same time, 
as we have seen above, the time interval between subse- 
quent reconnections becomes longer and longer. Eventu- 
ally, a point may be reached when the amount of plasma 
precipitated during 5t will become equal to the amount 
of material injected into the loop during an evaporation 
episode. After that point, there will be no further net 
secular change in the coronal density; the system will 
attain a stable limit cycle behavior. To evaluate the ex- 
act conditions characterizing this state one will need to 
calculate the amount of material evaporated following a 
given small flare and the rate at which the plasma returns 
to the photosphere as a result of gradual cooling. This 
would require a detailed model of the thermal structure 
along loop (e.g., along the lines of RTV78), but this can 
(and should!) definitely be done. For now, however, we 
just want to get some qualitative feeling. Therefore, we 
shall just say, tentatively, that the system enters this sta- 
ble cyclic regime (with no secular density gain) when the 
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radiative cooling time becomes equal to the time between 
reconnection events. This yields the following condition: 

Jt(n)=7-i :^=T,ad(n). (34) 

This equation can now be viewed as a condition that 
determines the long-term equilibrium plasma density in- 
side the loop, n*. Substituting equation (fT2)) for the 
radiative cooling time and equation (|27p for Be, we ob- 
tain as a function of the loop temperature Tioop, 7, L, 
and Bz'- 

^ 2 . 10^° cm-3 b'J! L''/' {jr^^^.o)'/' f ^) , 

V Q~22 J 

(35) 

where Q-21 = (9(Tloop)/(10 ^^ergcm^s ^) and where 
for convenience we defined Tiad.o as the radiative cool- 
ing time corresponding to nio — 1, Tioop = 1 MK, 
and (5-22 = 1: Trad,o — 400 sec. 

Correspondingly, the characteristic threshold recon- 
necting magnetic field is 

-So,* — Bc{n^,,L, Bz) 

45 G Ll" bT, (7r.ad,o)^/^ (^) , (36) 

and the characteristic time interval between subsequent 
reconnection events is 

At* =ri.ad(?i*) 

200 sec BJ,^/' lI^' (^r.ad.o)"^/^ (^) '1?7) 

The characteristic amount of energy released in each 
reconnection event comprising the limit cycle is 

B^ 

X (7r.ad,o)^/^ i^T'^^^^ 

and hence the time-averaged rate of magnetic dissipation 
in the loop is 

TP 

= — ^ ~ 4 • 10^6 erg sec-i x 

1 /3 

X LlgLg B^ 2 (TT^rad.o) 1 g I ,(39) 

Note that most of the above parameters have weak 
dependence on L, except for and H^, whose strong L- 
dependence is almost entirely due to the volume involved. 

To make any further progress, we need two more rela- 
tionships: (1) the cooling function Q{T) and (2) a scal- 
ing for the loop temperature Tioop. Note that both are 
subject to significant uncertainties. In particular, the ra- 
diative cooling function exhibits a complicated behavior 
in the relevant temperature range (1-10 MK); for exam- 
ple, according to Cook et al. (1989), it decreases sharply 
between T = 1 MK and T = 3 MK, and then stays 
nearly flat for T > 3 MK. Therefore, in principle, one 
should not expect simple power-law scalings to result at 
all. Furthermore, there are still significant disagreements 
in the literature regarding Q{T) (e.g., Raymond et al. 



1976; Raymond & Smith 1977; RTV78; Cook et al. 1989; 
Landi & Landini 1999; Aschwanden et al. 2000a). For 
these reasons, any realistic quantitative analysis of ra- 
diative cooling in the context of the present model is a 
highly non-trivial task, adding even more ambiguity to 
this already complicated picture. Getting into such a 
high level of sophistication and detail would exceed the 
overall level of accuracy of our model. We shall therefore 
leave it for a future study. 

In addition to the cooling function, we need an ex- 
pression for the loop temperature, Tioop, which strongly 
affects plasma cooling, in terms of the other loop parame- 
ters. In principle, the temperature, or, more generally, its 
distribution along the loop, is determined by the hydro- 
static balance in conjunction with the energy transport 
along the loop, which includes thermal conduction, ra- 
diative losses, and distributed heating. Such an analysis 
has been first performed by RTV78 for the simplest case 
of uniform and steady heating; as a result, they derived 
the famous RTV scaling for the loop-top temperature: 
Tj^TV = . io3(pL||)i/3^ where p is the loop pressure 
and T|| is its length (which is different and usually larger 
than our current-sheet length L; in our model, we shall 
regard it as constant). Now, it is true that there is a lot 
of disagreement in the modern literature regadring the 
applicability of the RTV scaling to the real solar corona 
and the observational support is doubtful (e.g.. Porter 
& Klimchuk 1995; Cargill et al. 1995; Kano & Tsuneta 
1995; Aschwanden et al. 2000b,2001a). This is not sur- 
prising taking into account that the RTV scaling was de- 
rived assuming stationary and uniform heating, and thus 
ignored the impulsive and intermittent nature of coronal 
energy dissipation (which is an intrinsic property of my 
model). In addition, to develop their theory, RTV78 re- 
lied on an old cooling function due to Raymond (see be- 
low), whereas somewhat different cooling functions have 
been used in recent years (Cook et al. 1989; Landi & 
Landini 1999; Aschwanden et al. 2000a). 

Nevertheless, despite its significant limitations, the 
RTV scaling has been highly infiuential in solar physics. 
It is still a highly-respected common standard, against 
which solar physicists often measure their theories and 
observations. Therefore, just to present an illustrative 
example, let us now combine the RTV scaling with our 
model. First, to be consistent, we need to adopt the 
cooling function due to Raymond (Raymond et al. 1976; 
Raymond & Smith 1977) that was used by RTV78:3 

qRTV(j,) ^ 2 • 10-22 T-'/' ergcm^ s'' , (40) 
apphcablc for 2-10 MK (see Appendix A of RTV78). 
Substituting this function into scalings ([55)1 - ([55]) . we get 

^ 10^° cm-3 L-^/^ tZ', (T^.ad.o)^/^ , (41) 

Bo, - 40 G b'J^,' lI" T^il, (7r.ad,o)^/^ , (42) 

AU ^ 160 sec B-^^' L^/^ T.Hl, hr,,^,,)-'/' , (43) 

El r in28 T r 2+2/9 nlO/9 rTilO/9 / \2//2a\ 

£;*~5 • 10 ''erg Tn^gig B^J T,Jpg (7Trad,o) (44) 
i?* ~ 3 • 10"^^ erg sec-i L^^^g bI^^^ x 

X ri^ofp,6 (7rrad,o)^/^ (45) 

^ Note that, to derive their scaUng, RTV78 actually used an 
approximation of the Raymond cooling function: Q{T) ~ T~^^^. 
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Next, recasting the RTV loop-temperature scaling in 
terms of the electron density instead of the pressure, 

T.^Jv^3.106Knl/^LjJ, (46) 

and substituting it into the above scalings, we get: 

T,«o"p^* - 7 . 10« K B^l^, L-'" l; J , (48) 
B^J-'^SOGB^'Lr'Ll^^j'Ji, (49) 
A^^^ ^ 900 sec Ll/'' lJJ jZl^' , (50) 

iJ^V ^ 2 5 . 1029 ^9/4 ^2+1/12 ^5/^3 ^3/^2 ^ ^^^^ 

xL;;rL^/^^<v^:i^ (52) 

where, in addition, we introduced 7-3 = 7 • 10^ sec and 
used Trad.o = 400 sec. 

A similar exercise could be conducted, for example, 
for an entirely different cooling function (Cook et al. 
1989; Landi & Landini 1999; Aschwanden et al. 2000a) 
and/or for a different scaling law [e.g., that derived by 
Cargill et al. (1995) or, observationally, by Kano & 
Tsuneta (1995)]. 

These estimates give us the characteristic scales of 
micro- and nanoflare-like energy-release events in a 
bright coronal loop of a fixed length fixed cross-loop 
spatial scale L, and a fixed axial magnetic field B^, sub- 
ject to continuous braiding at a fixed rate 7. It is of 
course understood that this is just an idealized model 
set-up and that, in practice, using such typical values 
straight-forwardly may not be meaningful because of the 
highly inhomogeneous nature of the solar corona. 

4. DISCUSSION: WHAT NEEDS TO BE DONE 

In this section I will discuss several open questions rel- 
evant to the physical picture presented in this paper — 
questions that I would like to see answered in the near 
future. Correspondingly, I will describe possible studies 
(mostly numerical) that I think ought to be performed 
in order to confirm, modiiy, or refute various elements of 
my model. 

First, the whole picture hinges, in part, on the premise 
that rcconncction is slow in the coUisional regime. Al- 
though in the above discussion I have for definiteness as- 
sumed that it is as slow as in the classical Sweet-Parker 
model, I don't actually think that this needs to be the 
case. If the true rate of classical resistive-MHD reconnec- 
tion turns out to be much faster than Sweet- Parker, the 
overall picture would still stand, qualitatively, as long 
this rate is still much smaller than the coUisionless re- 
connection rate. The only thing that matters is that 
there are two regimes, slow and fast, and the transition 
between them is governed by the coUisionality of the re- 
conncction layer. 

Nevertheless, from the point of view of my picture, 
studies of slow collisional (resistive MHD) reconnection 
are just as important as those of fast coUisionless recon- 
nection. One would like to tighten up the case for slow 
reconnection. Thus, I would like to encourage more stud- 
ies of resistive reconnection, to really confirm that it is 



slow. Thus, it is important to address the following spe- 
cific issues, best using numerical simulations: 

1) Does the presence of small-scale 3D MHD turbu- 
lence enhance reconnection rate, as suggested by Lazar- 
ian & Vishniac (1999)? If it does, how large is the en- 
hancement and what determines it? Could it be as fast 
as coUisionless reconnection, in which case it would con- 
stitute a challenge to the picture presented in this paper? 
To address this issue, one would have to perform a 3D 
MHD simulation involving a regular large-scale recon- 
nection layer with a super-imposed 3D MHD turbulence. 
Just as importantly, one needs to investigate the possibil- 
ity that 3D MHD turbulence is spontaneously generated 
inside the layer by the secondary-instability mechanism 
of Dahlburg et al., (1992, 2005), and evaluate its effect 
on the time-averaged reconnection rate. 

2) What is the effect of the 2D tearing instability inside 
the resistive Sweet-Parker layer, on the overall, time- 
averaged reconnection rate? It is expected that tearing 
instability may become important in reconnection layers 
that are very long, with aspect ratios L/5 > 100. It 
may then make the reconnection process inherently time- 
dependent, bursty. Therefore, in order to investigate this 
issue, one would have to perform only a 2D simulation, 
but with a very high resolution (corresponding to 5 > 
10^) and with a very long duration (to be able to average 
over many bursts). 

3) The majority of existing numerical studies of 
resistive-MHD reconnection use a constant uniform re- 
sistivity. It is important to check whether the main re- 
sults will be the same if one uses the actual Spitzer re- 
sistivity. This is especially important in the context of a 
magnetically-dominated environment, such as the solar 
corona, where the temperature inside the reconnection 
layer may be much higher than the ambient tempera- 
ture. As a result, the Spitzer resistivity may vary by a 
factor of a hundred across the layer, being much smaller 
at its center. 

4) One should further pursue numerical studies of mag- 
netic reconnection that include both resistive MHD and 
Hall eff'ect, along the lines of the recent work by Cas- 
sak et al. (2005, 2006, 2007). Of particular interest is the 
intermediate regime, when Sgp > di,pi but at the same 
time fie > Ve- The objective of such simulations would 
be to observe the transition from slow to fast reconnec- 
tion within one simulation run as the upstream plasma 
parameters, such as Bq orn, are gradually changed, pass- 
ing the critical point. Whereas the results published by 
Cassak et al. (2005, 2006, 2007) are already very impor- 
tant, the relevant parameter space is large and similar 
studies need to be extended to other regimes, most no- 
tably to a situation where /3upstroam ^ 1- Also, these 
results need to be confirmed by other groups, preferably 
with a more realistic electron-to-ion mass ratio, a larger 
box size, and higher numerical resolution. 

5) Future numerical studies of collisional reconnection 
should include realistic energy balance taking into ac- 
count both Ohmic heating and the electron thermal con- 
duction. The goal here would be to calculate the values 
of Tg and Ue at the center of the reconnection region. 
Also, one needs to study temperature equilibration be- 
tween ions and electrons. 

6) Effect of viscosity on resistive-MHD reconnection 
needs to be assessed. 
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7) The effect of a strong guide field on resistive-MHD 
reconnection and on the transition to the fast regime 
(e.g., Cassak et al. 2007) needs to be investigated in more 
detail. 

8) Finally, I would like to strongly encourage further 
experimental (laboratory) studies of coUisional reconnec- 
tion, especially in the large-S* limit. 

In addition to the above questions related to coUisional 
reconnection, there are several important issues related 
to collisionless reconnection: (1) What is its physical na- 
ture of anomalous resistivity due to wave-particle interac- 
tions and under what conditions is it excited? (2) What 
is the structure of the Petschek-like reconnection layer 
for a given functional form of anomalous resistivity and 
what is the resulting reconnection rate in terms of the 
basic physical parameters? Does the anomalous resistiv- 
ity spread along the separatrices (Petschek's shocks) or is 
it present only in the small central region? If the latter is 
the case, then how does the plasma crossing the shocks 
gets heated? (3) When does laminar a Hall-effect (or, 
in general, two-fluid) reconnection take place? How do 
Hall effect and anomalous resistivity co-exist and interact 
with each other? (4) How rapid is two-fluid reconnec- 
tion? What parameters affect the reconnection rate in 
this regime? (5) What is the effect of guide field on trig- 
gering and saturation of Hall-regime reconnection and 
on anomalous resistivity? (6) Is collisionless reconnec- 
tion process bursty and, if so, what is the time-averaged 
reconnection rate? (7) What is the overall partitioning of 
the released energy between the bulk kinetic energy, elec- 
tron and ion thermal energies, and nonthermal particle 
acceleration? 

Finally, in order to be realistic, future numerical simu- 
lations of the solar corona should include all of the follow- 
ing (c.f., Miyagoshi & Yokoyama 2003; Klimchuk 2006): 
(1) flux emergence processes and random motions of the 
field-line footpoints; (2) physically-motivated prescrip- 
tion for transition to fast reconnection, such as the one 
suggested in this paper; such a prescription would thus 
play a role of a subgrid model used in a large-scale MHD 
simulation of the corona; (3) mass exchange between the 
corona and the solar surface (e.g., chromospheric evap- 
oration and plasma precipitation); (4) optically-thin ra- 
diative energy losses and thermal conduction (including 
the contribution due to nonthermal electrons) along the 
magnetic field lines. 

5. CONCLUSIONS 

Magnetic reconnection research started 50 years ago 
in the field of Solar Physics, with the Sweet-Parker 
(Sweet 1958; Parker 1957) model for solar fiares, followed 
by the Petschek (1964) theory a few years later. These 
studies tackled the most fundamental micro-physical as- 
pects of the reconnection layer. Over time, however, the 
focus of solar reconnection research has shifted away from 
local basic physics of reconnection, despite the fact that 
reconnection has been confirmed observationally to be 
the key process responsible for magnetic energy release 
in flares (e.g., Tsuneta 1996, Yokoyama et al. 2001). 
Instead, the forefront of research on the fundamental 
physics of reconnection has moved in recent years to 
other fields, most notably, to Space Physics, where de- 
tailed in-situ measurements using Earth-orbiting space- 
craft are now available (e.g., Oieroset et al. 2001; Na- 



gai et al. 2001; Mozer et al. 2002), and to Laboratory 
Plasma Physics, where several dedicated experiments 
(Yamada et al. 1997; Egedal et al. 2007; Brown 1999) 
have already made fundamental contributions to our un- 
derstanding of reconnection. Most importantly, a lot of 
progress has been recently made using numerical simula- 
tions, again, mostly in the context of the Earth Magne- 
tosphere. 

The basic paradigm that emerges as a result of all 
these studies can be summarized as follows (see § 12. ip . 
The starting point is the realization that there are in- 
deed two reconnection regimes. The first one is a slow 
(Sweet-Parker) resistive-MHD regime that is realized in 
relatively dense, coUisional plasmas. The second one 
is a fast (Petschek-like) regime that takes place in col- 
lisionless plasmas. The mechanism for the fast colli- 
sionless reconnection can be either a locally-enhanced 
anomalous resistivity due to micro-instabilities triggered 
when the current density exceeds a certain threshold; 
or the Hall effect. In either case, one can formulate 
an approximate condition for the transition between the 
slow coUisional and the fast collisionless regimes. If 
the guide field is not much larger than the reconnect- 
ing magnetic field, this condition is 6sp < di, where 
6sp is the thickness of the Sweet-Parker reconnection 
layer and di is the collisionless ion skin depth. One 
can further rewrite this condition in terms of the clas- 
sical electron mean free path Ag^mfp inside the layer as 

where L is the global sys- 
tem size (Yamada et al. 2006). Due to the strong tem- 
perature dependence of Ae,mfp, this form of the condi- 
tion highlights the need to estimate the electron temper- 
ature Te at the center of the layer. Using considerations 
of pressure balance and energy conservation (see § 12. 3p , 
one can express Tg in terms of the reconnecting magnetic 
field i?o and the background plasma density n. Then, the 
collisionless reconnection condition can be recast in terms 
of L, Bq and n, see equation (fTOl) in § 12.21 In the case 
with a strong guide field i?^ ^ i?o, the corresponding 
condition is (5sp < pi and this leads to an additional fac- 
tor (Bq/Bz)'^ in the expression for the critical length Lc 
[see eq. (^51 ]. 

One of the main driving forces behind this paper is the 
author's desire to bring the recently-obtained knowledge 
about reconnection back to Solar Physics and to use it 
productively to build a better understanding of the solar 
corona. Although most of the present discussion is also 
relevant to solar flares, in this paper I focus predomi- 
nantly on the problem of solar coronal heating (see § [3]) . 

In the context of Parker's (1983, 1988) nanoflare the- 
ory of coronal heating, magnetic energy release in the 
solar corona takes place in the form of many unresolved 
reconnection events (nanoflares). One of the most im- 
portant features of this picture is the intermittent char- 
acter of energy release. Random footpoint motions lead 
to continuous twisting of elementary magnetic strands 
around each other, which, in turn, leads to the forma- 
tion of many small current sheets in the corona. Current 
sheets may form either in finite time, as suggested by 
Parker (1983, 1986), or exponentially in time, as was 
demonstrated by van Ballegooijcn (1986) and later nu- 
merically by Galsgaard & Nordlund (1996); for our pur- 
poses, it does not matter which one is correct. What 
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matters is that thin current layers do eventually form. 
Free magnetic energy accumulates for a while and is then 
suddenly released in distinct fast reconnection events. In 
the present paper, I suggest that the transition between 
the slow and fast reconnection regimes plays a key role 
in determining when a given nanoflare will take off and 
how much energy will be released (see also Cassak et al. 
2005). My model can thus be regarded as an alterna- 
tive to the secondary-instability mechanism proposed by 
Dahlburg et al. (2003). Furthermore, I argue that the 
fact that the fast reconnection condition involves the am- 
bient plasma density is an important part of the story. 
The reason for this is that the density in the corona is not 
a fixed constant; in a given loop, it constantly changes in 
response to radiative cooling and chromospheric evapo- 
ration caused by coronal energy-release events. The ba- 
sic picture here is the following. A coronal energy re- 
lease leads to an increase in density, thus making the 
plasma more collisional. This temporarily inhibits fast 
reconnection in the given region until the density de- 
creases again (on the radiative cooling timescale) to be- 
low a certain critical value. At this point fast reconnec- 
tion again becomes possible. On the longer time-scale, 
a quasi-periodic behavior is established, characterized by 
repeated cycles that include fast reconnection events, fol- 
lowed by chromospheric evaporation episodes, followed 
by relatively long (^ 1 hour) and steady periods during 
which free magnetic energy in the loop builds up and the 
plasma gradually cools down. Thus, coronal heating can 
be viewed as a self-regulating process that statistically 
keeps the density roughly near the critical value for the 
fast reconnection transition. In other words, the system 
constantly fluctuates around the state of marginal coUi- 
sionality as defined by the coUisionless reconnection con- 
dition. In the long-term statistical equilibrium, a balance 
is maintained in which the amount of plasma pumped 
into a coronal loop as a result of an evaporation episode 
is equal to the amount of plasma drained down to the sur- 
face during the gradual radiative cooling stage that takes 
place between two subsequent fast-reconnection events. 
The characteristic equilibrium density, the time interval 
between reconnection episodes, and the energy released 
in each such episode can be estimated in terms of the 
loop's longitudinal magnetic field, its characteristic size, 
and the footpoint driving rate (see § ^ . 

Finally, I believe that the physical framework de- 
veloped in this paper should also be applicable to 
magnetically-dominated coronae of other astrophysical 
objects, such as other stars and accretion disks (Good- 
man & Uzdensky 2007). 
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tions. 
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